function [ER_Mean] = ER_Mean_Prism(Data_Summary, animalName, drug)

for n = 1:size(Data_Summary, 2)
    if strcmp(Data_Summary(n).animalName, animalName) && strcmp(Data_Summary(n).drug, drug)
        ER(n,1) = Data_Summary(n).Session1_Mean_Ca_Event_Rate;
    else
        ER(n,1) = nan;
    end
end
ER_Mean = mean(ER, 'omitnan');

end